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Abstract 



A new method has been developed to optimize, in terms of aerodynamic wave drag 
minimization, arbitrary (nonaxisymmetric) hypersonic vehicles in modified Newtonian flow, while 
maintaining the initial volume and length of the vehicle. This new method utilized either a surface 
fitted Fourier series to represent the vehicle's geometry or an independent point-motion algorithm. 
In either case, the coefficients of the Fourier series or the spatial locations of the points defining 
each cross section were varied and a numerical optimization algorithm based on a quasi-Newton 
gradient search concept was used to determine the new optimal configuration. Results indicate a 
significant decrease in aerodynamic wave drag for simple and complex geometries at relatively low 
CPU costs. In the case of a cone, the results agreed well with known analytical optimum ogive 
shapes. The procedure is capable of accepting more complex flow field analysis codes. 


Nomenclature 




X 


y> z 

s 


s 

Y 

A 


= surface pressure coefficient 
= stagnation pressure coefficient 
= static pressure at a point 
= free stream static pressure 
= free stream Mach number 

= angle between free stream and normal to the surface of a vehicle 
= coefficients of Fourier trigonometric series for coordinates at cross section i 
= Cartesian coordinate along the axis of the body 
= Cartesian coordinates of a contour point at cross section i 
= body cross section contour-following coordinate 

= normalized body cross section contour-following coordinate 

= specific heat ratio of the gas 
= area of a panel on the body surface 


F = aerodynamic force applied to a panel 

n = unit normal to the body surface 

M = number of terms in the Fourier trigonometric series 

S = least squares summation 

FAC = percentage change in design variable 


Subscripts 
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i = i f h cross section of the vehicle 

j = j 1 * 1 point of a cross section contour 

m = m^ 1 coefficient of a Fourier trigonometric series 

n = angle between free stream and local body surface normal 

00 = free stream value 

Introduction 

Although optimization of axisymmetric hypersonic bodies has been accomplished in the 
past [1,2], the aerodynamic drag minimization of an arbitrary hypersonic vehicle has not been 
attempted [3]. The objective of this paper is to present an optimization procedure for arbitrarily 
shaped hypersonic vehicles. While there are certainly some limitations in this paper, including the 
choice of flow field solver and non-convergence for some shapes, it demonstrates that optimization 
of numerous variables can indeed be done and that this can be applied to complex configurations. 

In hypersonic flow (M„ > 5.0), the flow around an object may be modeled using an impact 
theory. In this theory, oncoming particles strike, or impact, the surface of the object and impart the 
normal component of their momentum to that body. Classical Newtonian theory has been shown 
to approach reality when the free stream Mach number approaches infinity and the value of the ratio 

of specific heats approaches y =1 [4], In the case of lower hypersonic Mach numbers, modified 
Newtonian theory has been shown to be quite satisfactory for predicting the aerodynamic forces 
and moments on a body [5]. Modified Newtonian theory has the main advantage of being 
extremely simple, accurate [6], and fast when faced with the thousands of flow field calculations 
needed in an optimization problem of this scope. Because of the use of modified Newtonian 
theory, it was implicitly assumed that the flow field was inviscid. 

In this study, modified Newtonian impact flow theory was used with a modified 
Newtonian constrained search optimization routine [7] to obtain vehicle shapes which had 
significantly lower wave drag in inviscid hypersonic flow. In the first part of the study, cross 
section coordinates of the body were represented with curve-fitted Fourier series. Curve-fitted 
Chebyshev series [8] were initially considered, but it was found that the Fourier series represented 
complex shapes, such as a Space Shuttle configuration, better than the Chebyshev series. The 
coefficients of the Fourier series, one set representing the y coordinate and one set for the z 
coordinate (Fig. 1), then became the design variables that were fed to the optimization routine. The 
optimization routine sequentially perturbed each of the coefficients by a small amount and 
determined the new shape that reduced wave drag while keeping the volume and length of the 
vehicle constant. In the second part of the study, the y and z coordinates of the vehicle s cross 
section were used as the design variables directly. Again, the optimization routine perturbed 
separately each of the coordinates at each of the cross section contour points. Then, it combined 
the changes into a new shape with lower wave drag while still honoring the constraints of constant 
volume and constant length of the vehicle. 

Numerical Models 

The first part of this investigation used a least squares Fourier series curve fit to represent 
the y and z coordinates of each half cross section, i , that is, 

M 

yjj = X A rni COS[(m-l)7CS' i j] 


(1) 
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M 

z u = X B mi sin[m7ts' ij ] 


m « 1 


where Sjj is a normalized contour-following coordinate (Fig. 1) such that 
s i,j = S i,j-1 + V (yi,j-yi.j-i ) 2 + (Zij-Zi.j-i ) 2 and s 10 = 


J i,jmax 


and A^j is determined from the least squares fit of the Fourier series 


Si = 


jmax 

X M 

■i 

m = 1 

J = 1 


dS, 


(A mi cos[(m-l)7ts ltJ ]) - yij r and = 0 


( 2 ) 


(3) 


(4) 


The coefficients B mi are determined in a similar way. Since it was assumed that the vehicle had a 
vertical plane of symmetry, the z-coordinates of the first (j=l) and the last (j=jmax) point of each 
half cross section were always zero thus ensuring symmetry across the y-axis. 

The local surface pressure coefficient, C p ;j , was calculated by the use of modified 
Newtonian impact flow theory, which states that 

Cp.ij = Cpo COS 0 n jj (5) 

where 0 n ^ is the angle between the free stream and the normal to the surface. The stagnation 
pressure coefficient, C^, is given by 


C = 

^po 


[( 


y + 1 


YM^ 2YM,,, - y + 1 


t 

-)y-1 


Y+l 


M 


JL 

2 >' i]. 


( 6 ) 


The pressure on a given segment of the body may then be calculated from the rearranged formula 
for Cp y, that is 


Pij P<» + 2 ^"P-ij y P°° M ~ 

The aerodynamic force on each surface panel is found by 
B ij Pij A jj tljj . 

so that the resultant force is obtained by summing up all of the panel forces 

B total = ^Lf B ij 

Aerodynamic wave drag was then the x-component of the resultant aerodynamic force. 


(7) 


( 8 ) 


( 9 ) 
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The optimization algorithm perturbed each of the Fourier coefficients to obtain a sli°htlv 
different shape. After pert rbing all of the coefficients and analyzing these new perturbed shapes, 
the optimization algorithm combined the changes into a new shape that met the constraints of 
constant volume and constant length, but which had a reduced aerodynamic wave drag. 

The second part of the study was exactly the same as the first except that instead of 
working with the Fourier series coefficients we worked with the v and z coordinates of the cross 
sections’ points directly. 

Results 


Four test cases were run for each part of this study. They consisted of a straight cone 
having circular cross section shapes, a straight cone having a four pointed star as a cross section, a 
stubby-wing shaped body and a Space Shuttle-like configuration. All cases were run at an angle of 

attack of 0°, a specific heat ratio of y= 1.4 and a free stream Mach number of M oc = 10. Notice that 

the values for y and M„ appear in Cpo which may be factored out of the pressure coefficient ratio. 
They affect the numerical amount of drag, but not the qualitative amount of drag. The x-axis for 
each case was chosen to coincide with the long axis of the body. The y and z-axes were then 
mutually perpendicular to the x-axis. 


Fourier Series Algorithm; For the initial part of the study, twenty terms in a Fourier 
senes for y and z coordinates were used for seven cross sections. Twenty terms were chosen 
because of the constraints of computational facilities (an IBM 3090 was used) and because twenty 
terms were able to represent the geometries of all four test cases, including the complex Space 
Shuttle shape. Only six of the cross sections were allowed to deform; the nose cross section was 
kept constant to serve as a tip. Thus there were 6 x 20 x 2 = 240 design variables. Twentv-five 
points were used per cross section; thus, the half body was discretized into 6 x 24 = 144 panels. 
FAC, the percent perturbation of A mi and B mi in the optimization algortihm, was set equal to 5%. 

For the case of a right circular cone (Fig. 2), after a total of 43 iterations, the program 
converged to an ogive configuration that had 47.96% less wave drag than the original conical 
configuration. Note that horizontal and vertical symmetry were maintained. 

The next shape tested was a four pointed star configuration (Fig. 3). The aerodynamic 
wave drag of this shape was reduced by 39.16% after 43 iterations. This case did not converge 
due to tishtailing of the fins and was terminated just before such fishtailing occurred. Note the 
streaking near the nose and the thinning of the points on the last cross section. Also, notice that 
vertical and horizontal symmetry was maintained and the fins gave an ogiving contour 

The third shape optimized was that of a "stubby wing" configuration (Fig. 4). The wave 
drag was reduced by 64.62% when the algorithm converged after 53 iterations while preserving 
the cross-axis symmetry. Again, note the streaking toward the nose and the smooth appearance of 
tour small fins along the wing tip line. 

The fourth test case was that of a Space Shuttle configuration (Fig. 5). After 22 iterations 
the W1 "§ su rf a ce crossed itself and the process was terminated. Aerodynamic wave drag was 
reduced by 18.52%. With careful scrutiny, one can notice that the centerbody has become ogived 
the wing thickness has been reduced, the wing roots have become filleted, and the underside of the 
fuselage has been reduced in size. 

The convergence histories (Fig. 6) indicate that the general trend is a monotonic decrease in 
drag, This trend can also be seen from the drag plot (Fig. 7) which shows the percentage of 
original drag remaining at a given iteration number (Table 1 ). 


EPIPt-MfltiQn Algorithm : For the point-motion algorithm, 21 points per half cross 

section were used. Only six cross sections were analyzed due to computer storage limitations 

^“ h inn y !l Ve Cr ° SS f Ct °u S , b r dng 3Ctive ’ this yielded 2 x 5 x 21 = 210 design variables and 5 x 
20 100 surface panels per half cross section. FAC for this algorithm was set to 0. 1 %. 
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Numerical optimization of a straight circular cone resulted in an ogive shape with 
axisymmetry successfully maintained. After 49 iterations, the shape converged to that of Fig. 8 
with a corresponding decrease in drag of 45.37%. 

The second shape was the four pointed star. With the point-motion algorithm, the shape 
converged after 35 iterations to that pictured in Fig. 8. There was a reduction of drag of 34.65%. 
As in the case of the circular cone, symmetry was maintained across both the y and z-axes. The 
resulting shape is very similar to that obtained by the coefficient algorithm. Considerably less 
streaking can be seen near the nose of the star, while the fin planforms exibited significant ogiving. 

The next shape optimized was the "stubby wing". The fin tips, after 24 iterations, crossed 
themselves and the process was terminated. However, a decrease of 40.58% of the original wave 
drag was achieved just before the shape cross-over. Note the development of the fins along the 
side of the vehicle (Fig. 8) and streaks near the nose, somewhat similar to those developed when 
using the coefficient algorithm. Once again, symmetry was maintained. 

The final case for the point motion algorithm was the Space Shuttle vehicle (Fig. 8). 
Similar to the case in the coefficient algorithm, the wing eventually crossed itself and the run was 
terminated. A decrease in drag of 27.38% was found after 16 iterations preceding the cross-over. 

For the point-motion algorithm, the convergence histories (Fig. 9) indicate monotonic 
decrease in wave drag for all four test configurations. Figure 10 demonstrates the total reduction in 
wave drag for the four test configurations when using the point-motion algorithm (Table 2). 

A comparison (Figs. 11 and 12) of the numerically optimized ogive shapes with 
analytically optimal ogives obtained by Sears and Haack and by von Karman [9] demonstrates the 
reliability and accuracy of the numerical optimization algorithms. 

Conclusions 

Two procedures, a coefficient algorithm and a point-motion algorithm, for aerodynamically 
optimizing aribitrarily shaped hypersonic vehicles have been shown to significantly reduce 
aerodynamic wave drag while keeping the vehicle's volume and length constant. Both 
formulations are very fast only because a modified Newtonian flow theory was used as the flow 
field analysis algorithm. These formulations would be very effective as preliminary design tools 
for unconventional hypersonic vehicles. The point motion algorithm can be used to keep parts of 
the original vehicle fixed, such as cabin size or wing thickness, during the optimization. More 
sophisticated flow field solvers that include viscosity and the effects of heat transfer could be 
substituted in place of the modified Newtonian theory during the final stages of the optimization. 
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Coefficient Test Case 

Drag Reduction (%) 

# of Optimization 
Cycles 

# of Analysis Calls 

Cone 

47.96 

43 

10493 

Star 

39.16 

43 

10493 

Stubby wing 

64.62 

53 

12689 

Space Shuttle 

\ 18.52 

22 

3T25 ( 


Table 1. Drag reduction, number of optimization cycles and analysis calls for the coefficient test 
cases that uses Fourier series representation. 
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a) initial 


b) final 



Figure 2. Straight circular cone shape; coefficient algorithm: a) initial shape, b) final shape. 



a) initial 


b) final 


Figure 3. Four pointed star shape; coefficient algorithm: a) initial shape, b) final shape. 



a) initial b) final 


Figure 4. Stubby wing shape; coefficient algorithm: a) initial shape, b) final shape. 
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Figure 6. Convergence histories for four test shapes when using coefficient algorithm; 
percentage change in wave drag per iteration. 



Figure 7. Convergence histories for four test shapes when using coefficient algorithm; 
total change in wave drag per iteration. 




Figure 8. Optimized shapes obtained with the point-motion algorithm. 
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Figure 9. Convergence histories for four test shapes when using point-motion algorithm; 
percentage change in wave drag per iteration. 



Figure 10. Convergence histories for four test shapes when using point-motion algorithm; 
total change in wave drag per iteration. 
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Ogives 



Figure 1 1 . Comparison of analytically optimized ogives and a numerically optimized ogive using 
coefficient algorithm. 


Ogives 



Figure 12. Comparison of analytically optimized ogives and a numerically optimized ogive using 
point-motion algorithm. 





